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ABSTRACT 

We present a large sample of 803 star-forming luminous compact galaxies 
(LCGs) in the redshift range z = 0.02 - 0.63, selected from Data Release 7 of 
the Sloan Digital Sky Survey (SDSS). The global properties of these galaxies 
are similar to those of the so-called "green pea" star-forming galaxies, in the 
redshift range z = 0.112 - 0.360 and selected from the SDSS on the basis of their 
green color and compact structure. In contrast to green pea galaxies, our LCGs 
are selected on the basis of both their spectroscopic and photometric properties, 
resulting in a ~ 10 times larger sample, with galaxies spanning a redshift range 
> 2 times larger. We find that the oxygen abundances and the heavy element 
abundance ratios in LCGs do not differ from those of nearby low-metallicity blue 
compact dwarf (BCD) galaxies. The median stellar mass of LCGs is ~ 10^ Mq. 
However, for galaxies with high EW(H/3), > lOOA, it is only 7x10^ Mq. The 
star formation rate in LCGs varies in the large range of 0.7 - 60 Mq yr~^, with 
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a median value of ~ 4 Mq yr~^, a factor of ~ 3 lower than in high-redshift 
star-forming galaxies at z > 3. The specific star formation rates in LCGs are 
extremely high and vary in the range ~ 10~^ - 10~^ yr~^, comparable to those 
derived in high-redshift galaxies. 

Subject headings: galaxies: abundances — galaxies: irregular — galaxies: ISM 
— galaxies: star formation — H II regions 



1. Introduction 



Recently, iCardamone et al. I (|2009[ ) reported on a new class of luminous compact galax- 
ies at redshifts z = 0.112 - 0.360, discovered by the Galaxy Zoo project from the Sloan 
Digital Sky Survey (SDSS) images. The Galaxy Zoo project collected simple morphological 
classifications of nearly 9 00 000 galaxies drawn fro m the SDSS, contributed by hundreds of 
thousands of volunteers (ILintott et al. II2008I . l201ll ). In particular, the attention was drawn 
on compact objects named "green pea" galaxies. Their appearance is mainly caused by the 
very strong [O III] A5007A optical emission line, with an unusually large equivalent width 
of up to ~ 1000 A, redshifting into the SDSS r band, resulting in a green color in the gri 
composite S PSS images (green being the color represented by the r band in these compos- 



ite images). ICardamone et al. I (120091 ) discussed a well-defined sample of 251 color-selected 
objects located mainly in low-density regions. They found that the median environmental 
density around the green pea galaxies is less than two-thirds of that around normal galaxies. 
Most of the green pea galaxies are strongly star-forming, although there are some active 
galactic nuclei interlopers, including eight newly discove r ed na rrow-line Seyfert 1 galaxies. 



For a further spectroscopic analysis, ICardamone et al. I (120091 ) selected a subsample of 80 
star-forming galaxies and found that they are low-mass galaxies (M ~ 10^'^ - 10^° Mq), with 
high star forma t ion ra tes (~ 10 Mq yr~^) and solar metall icities, if the solar c alibra tion of 



Asplund et al. I (120091 ). 12 -|- log (O/H)© = 8.7, is adopted. ICardamone et al. I ((20091) found 



that green pea galaxies have some of the highest specific star formation rates seen in the local 
Universe (up to 10~^ Mq yr~^), and are similar in size, mass, luminosity and metallicity to 
luminous blue compact galaxies. They are also similar to high-redshift ultraviolet-luminous 
galaxies (Lyman-break galaxies and Ly-a emitters). Those authors suggested that the green 
pea galaxies may be the last remnants of a mode of star formation common in the early 
Universe, and therefore provide an excellent local laboratory for studying the extreme star 
formation processes that occur in high-redshift galaxies. 



Later, lAmorm et al. I (l2010l) investigated the oxyge n and nitrogen chemical abundances 



in the green pea galaxies of ICardamone et al. I (120091 ). They found these systems to be 
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genuine metal-poor galaxies, with mean oxygen abundances 12 +log 0/H ~8.0, or 0.2 solar. 
While the N/0 r atios of green pea galax ies follow the relation with stellar mass of local star- 
forming galaxies, lAmorm et al. I (|2010( ) find that the green pea mass-metallicity relation is 
offset by ~ 0.3 dex to lower metallicities. They argue that recent interaction-induced inflow 
of gas, possibly coupled with a selective metal-rich gas loss, driven by supernova winds, may 
explain their findings. Interaction-induced infall of gas would also be consistent with the 
known galaxy properties such as high specific star forrn ation rates, extreme compactness, and 
disturbed optical morphologies. I Amor in et al. I (120101 ) concluded that the green pea galaxy 
properties seem to be not common in the nearby universe, suggesti ng that the green pea phase 



is a short and extreme stage of their evolution. They agreed with ICardamone et al. I (|2009[ ) 



that these galaxies may allow to study in great detail many processes, such as starburst 
activity and chemical enrichment, under physical conditions approaching those in galaxies 
at higher redshifts. 

In this paper, we wish to further th e stud ies of green pea galax i es an d reexamine the 
above conclusions of ICardamone et al. I ( 120091 ) and lAmorm et al. I (120101 ) by studying a 
considerably larger sample of luminous compact emiss ion-line galaxies selecte d from the 
Data Release 7 (DR7) of the Sloan Digital Sky Survey ( lAbazajian et al. II2009I ). We adopt 
a n entirely different approa ch in the selection of our galaxy sample. In contrast to the study 
of ICardamone et al. I (|2009[ ). our luminous compact galaxies are selected using spectroscopic 
criteria which are supplemented by morphological criteria. The selection criteria of our 
galaxies are described in Section 2. Their element abundances are derived in Section 3. 
The luminosity-metallicity relation for the sample galaxies is discussed in Section 4. In 
Section 5 we describe a technique for galaxy stellar mass determination, taking into account 
the gaseous continuum emission. We derive stellar masses for the young and old stellar 
populations in each galaxy. In Section 6 we compare the physical properties of starbursts in 
luminous compact galaxies. The mass-metallicity relation for the sample galaxies is derived 
in Section 7. In Sections 8 and 9, we examine respectively their star formation rates and the 
age of their stellar populations. Our conclusions are presented in Section 10. 



2. A sample of luminous compact galaxies 
2.1. Selection criteria 



The photometric selection criteria used by ICardamone et al. I (120091 ). based on the green 
appearance of a galaxy in the SDSS images, pick out galaxies in a relatively narrow range of 
redshifts, z ~ 0.1 - 0.3: those redshifts are such as to make the strong H/3, [O ill] A4959A 
and [O III] A5007A optical emission lines fall in the SDSS r band. Such an approach fails to 
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find galaxies with similar properties at lower and higher redshifts because their colors would 
then not be green. 

To avoid this restriction, we adopt here a different approach. We select luminous com- 
pact emission-line galaxies (LCGs) by using criteria that are based on both spectra and 
images in the SDSS DR7. These criteria are as follows: 

1) the extinction-corrected luminosity of the H/3 emission line is greater than L(H/3) = 
3x10^° erg s~^. The extinction coefficient is derived from the observed Balmer decrement in 
the SDSS spectra. 

2) the equivalent width of the H/3 emission line is high, EW(H/3) > 50 A. This criterion 
helps to select only objects with strong emission lines in their spectra (Fig. [T]), and thus 
cont aining young starbu r sts wi th ages 3-5 Myr. We note that most of the green pea galaxies 



from ICardamone et al. I (120091 ) have EW(H/3) > 100 A. Thus, our sample of LCGs covers a 



larger range of EW(H/9). 



Contrary to ICardamone et al. I ( l2009l ). we do not use the equivalent width of the [O ill] 
A5007 emission line as a selection criterion. This is because the strength of the [O ill] A5007 
emission line in young starbursts is metallicity-dependent. It is strongest at 12+logO/H ~ 
8.0, decreasing at both higher and lower O abundances. Therefore, using the strength the [O 
III] A5007 emission line as a selection criterion would unduly bias the sample toward galaxies 
with a narrow range of O abundances near 12-|-logO/H ~ 8.0. On the other hand, the use 
of the H/3 emission line allows us to pick out galaxies with a considerably larger range of O 
abundances. 

3) only galaxies with a well-detected [O III] A4363 A emission line in their spectra, with 
a flux error less than 50% of the line flux, are selected. This criterion allows an accurate 



abundance determination using the direct method, in contrast to ICardamone et al. I (120091 ) 
who used the less accurate empirical strong-line method for abundance determination. If we 
do not impose criterion 3, but only criteria 1 and 2, then the number of selected galaxies 
would increase only by ~20%. Thus, when criterion 3 is imposed, the decrease in the number 
of galaxies is small enough so as not to introduce large biases in the physical properties of 
the selected sample. On the other hand, the benefit of being able to determine accurate O 
abundances and study trends with metallicity is inestimable. 

4) galaxies with obvious evidence of Seyfert 2 spectral features, such as strong [Ne v] 
A3426, He ll A4686, [O l] A6300, and [S ll] A6717, A6731 emission lines are excluded. Thus, 
we select only star-forming galaxies. 



5) on their SDSS images, galaxies are nearly compact at low redshifts, and unresolved 
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at high redshifts. 



Using these criteria, we constructed 4 subsamples of LCGs: 1) subsample 1 includes 
161 galaxies with EW(H/9) > lOOA and a round shape, without evident signs of disturbed 
morphology. Thus, this su bsample should be most s imilar in properties to the star- forming 
green pea galaxy sample of lCardamone et al. I (120091 ): 2) subsample 2 consists of 159 galaxies 
with EW(H/3) > lOOA and with some sign of disturbed morphology suggesting interaction 
and the presence of more than one star-forming region; 3) subsample 3 is the same as 
subsample 1, but it consists of 116 galaxies with a lower EW(H/9), covering the range lOOA 
> EW(H/3) > 50A; and 4) subsample 4 is the same as subsample 2, but it consists of 367 
galaxies with a lower EW(H/3), lOOA > EW(H/3) > 50A. The galaxies from subsample 
4 often show two or more compact knots of star formation. Thus, the total number of 
selected ga laxies is 803, or ^ 10 ti mes larger than the sample of star-forming green pea 
galaxies of lCardamone et al. I (120091 ). In general, the selected LCGs are faint objects with 
low extinction. Their median SDSS g apparent magnitude and reddening E{B — V) range 
from 20.3 and 0.10 for subsample 1 to 18.4 and 0.15 for subsample 4 (Tabled]). 



2.2. A Luminosity Upper Limit 

The distribution of the H/? luminosity L{Y{(5) with redshift z is shown in Fig. HJ It is seen 
that the selected galaxies span a large range of redshifts, from ~ 0.02 to ~ 0.63. The apparent 
increase of L(H/3) with redshift, for 2;<0.2, is likely the result of a selection effect (called the 
Malmquist bias): only the brightest galaxies at high redshift were observed spectroscopically 
in the SDSS. Fig. Hjshows that galaxies with EW(H/3) > lOOA are distributed more uniformly 
with redshift than those with lower EW(H/3), which tend to populate the lower redshift 
region. This difference is likely also due to a selection effect. The [O ill] A4363A emission 
line in objects with a high EW(H/3) is stronger. Hence it can be measured with a good 
accuracy even in high redshift objects, which helps to fulfill the selection criterion 3 above 
and allows the inclusion of the galaxy in the LCG sample. 

Fig. |2] shows a remarkable feature: there is an upper limit of L(H/3) at the level of 
~ 2.5x10^^ erg s~^ (shown by the horizontal dashed line), which is almost not redshift 
dependent, despite the fact that none of our selection criteria sets an upper limit to the 
H/3 luminosity. The presence of this upper limit suggests the existence of a self-regulating 
mechanism in star formation which somehow forbids the formation of more luminous star- 
forming regions. There is probably a negative feedback mechanism due to the intense UV 
radiation of the present starburst which prevents further star formation in LCGs. While 
beyond the scope of the present paper, the role of the ionizing radiation from young massive 
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starbursts in shutting off star formation beyond a given tfiresliold luminosity, sliould be 
investigated. 



2.3. Blue compact dwarf, green pea and luminous compact galcLxies: a 

continuum in luminosity 

The range of ^ in Fi g. [2] is ~ 2.5 times la rger than that f or the green pea galaxies of 



Cardamone et al. I (|2009[ ). Those authors and lAmorm et al. I ( 120101 ) suggested that green 



pea galaxies form a new class of objects. We will argue now that this is not the case and that 
they are just a subset of LCGs. Nearby {z ~0) luminous compact emission-line galaxies have 
been known for decades. Examples of such objects are, in order of decreasing L(B.(3) (we 
give the logarithm of the aperture- and extinction-corrected L{}if3) in parentheses for each 
object in units of erg s-^): Mrk 930 (41.57), SBS 0335-052E (41.24), HS 0837+4717 (41.10) 
and II Zw 40 (40.81). Some of them have been studied in great detail (see e.g. the studie s 
of SBS 0335-052E by llzotov et al. lIlQQoi Il997l : iThuan et al.1ll997t iThuan fc Izotovll2005h . 
These objects are known as blue compact dwarf galaxies (BCD) a nd are plotted as large 
open circles in Fig. |2l A recent review of BCD properties is given bv lThuanl f l2008h . We see 
that the -L(H/3) of these BCDs are comparable to those of LCGs at the low-luminosity end 
of the sample. Other well-known BCDs have lower L(H/3), below our cut-off: Toll214-277 
(40.43), I Zw 18 (39.72 for a distance of 15 Mpc) and Mrk 209 (39.55). In other words, LCGs 
are the same type of objects as BCDs except that they are more luminous in the mean. 

By the same token, green pea galaxies do not form a distinct class of objects from 
LCGs. They are just a subsample of the class of LCGs, in a restricted redshift range. At 
other redshifts, the same galaxies would not look green, but blue, pink or brown. This is 
illustrated in Fig. |3] where we show SDSS images of eight galaxies from our LCG sample, 
located at different distances. LCGs in the upper panel are galaxies from subsample 1 with 
round shape, while those in the lower panel are galaxies from subsample 2 with extended 
emission and/or disturbed morphology. It is seen that, as the strong emission lines (the 
Balmer lines and the [O ill] A4959A, A5007A lines) shift through the different SDSS filters, 
the colors of the LCGs as seen in the composite SDSS gri images change with redshift. They 
are blue at low redshifts, becoming pink at 2; < 0.1, green at 2 ~ 0.2 and brown at ^ ~ 0.6. 

We next compare the colors of the LCGs to those of normal galaxies and QSOs. We 
show in Fig. H] the g — r vs r — i color-color diagram for LCGs (filled and open circles and 
triangles). For comparison are also shown representat ive objects classified in the SDSS as 



galaxies (black dots) and QSOs (ISchneider et al. Il2010[ ) (grey dots). The g~r and r — i colors 



for LCGs span very large ranges, extending over 2.5 mag and 3.5 mag respectively. These 
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ranges are the largest for LCGs with a high H/3 equivalent width EW(H/3) > lOOA (filled 
symbols). Such large color ranges are not due to differences in galaxy physical properties, 
but to a large redshift range. Galaxies with strong emission lines at z ~ 0.1 - 0.3 would 
appear in the left upper corner of the diagram, while the same galaxies at z ~ 0.4 - 0.5 
would occupy the region with r — i > 1.0, in the extreme right part of the diagram. The 
location of the LCGs in the g — r vs r — i diagram is very different from that of the bulk 
of galaxies without or with weak line emission in their spectra (grey dots). It is interesting 
to note that many LCGs are located in the region with {g — r, r — i) ~ (0.0 - 0.5, 0.0 - 
0.5) populated by QSOs. A significant fraction of these LCGs are automatically classified 
in the SDSS as QSOs because of their compact morphology, while their spectra are that of 
star-forming galaxies rather than AGNs. 



What is the o verlap between the 



green 



pea sample of ICardamone et al. I (|2009[ ) and 



our LCG sample? ICardamone et al. I ( l2009l ) selected star-forming green pea galaxies by 
requiring them to be in the left upper region of the {g — r, r — i) diagram, delineated by the 
two straight lines in Fig. [H We have encirc l ed in this figure all the 66 LCGs that are also in 
the green pea sample of ICardamone et al. I ( l2009l ). The latter sample contains a total of 80 
objects. Out of the 14 missing objects, 12 were rejected because of their noisy spectra with 
a non-detected or weak [O ill] A4363 emission line, although they were present in the initial 
sample. Thus, our LCG selection criteria missed only two green pea star-forming galaxies 
out of a total of 80 galaxies, i.e. the overlap factor between the green pea and LCG samples is 
98%! On the other hand, more than half of the LCGs in the regio n of green pea star-formiri g 



galaxies delineated by the two straight lines were not selected by ICardamone et al. I ( 120091 ). 
All these galaxies have a green color on the SDSS images and are at redshifts z ~ 0.1 - 0.3. 
As for the remaining LCGs, to the right of the solid lines in the color-color diagram, they 
have redshifts z < 0.1 or z > 0.3, and would not be picked up by the green pea selection 
criteria, although they are the same type of objects, but at different redshifts. 



2.4. Low-metallicity AGNs in Green Pea Galaxies 



It is interesting to note that the H/5 luminosities of the five known low-metallicity 
compact galaxies that are thought to contain AGNs, as evidenced by the strong broad 
components of their Ha and H/3 emission lines (llzotov fc Thuan Il2008t llzotov et al. Il2010l ) 
put them in the class of LCGs. These objects fulfill all the selection criteria discussed above 
for LCGs. They are shown by the large filled circles in Fig. |2j Their redshifts are in the 
range 0.1-0.3, and therefore they can also be classified as green pea galaxies. However, in 
the {g — r) - (r — i) color-color diagram in Fig. |U only two low-metallicity AGNs fall into the 
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green pea region. The two other low-metalhcity AGNs are located outside both the green 
pea and LCG regions. The latter two objects show abnormally red g — r an d r — i colors. 



proba bly due to high dust extinction. The fifth low-metallicity AGN studied by llzotov et al. 



( 120101 ) is not shown in Fig. H] because it is beyond the SDSS sky area and therefore we have 



no g,r,i colors for it. 

It is not unreasonable that AGNs tend to occur in more luminous and hence more 
massive compact galaxies, as these have more massive star formation and hence are more 
favorable to the formation of intermediate-mass black holes. It would be interesting to obtain 
high-resolution and high signal-to-noise ratio spectra of more objects in the LCG sample to 
look for more objects with broad components which may contain AGNs. 



2.5. The [O ill] A5007/H/3 vs. [N ll] A6583/Ha diagram 



Figure O shows by small filled and open symbols the position s of LCGs in the [O 



III] A5007/H/3 vs. [N ll] A6583/Ha diagram f lBaldwin et al. I Il98ll . hereafter BPT). The 



star and filled square represent respectively the well-studied BCDs SBS 0335-052E and HS 
0837+4717. HS 0837+4717 lies in the same region as the LCGs while SBS 0335-052E lies to 
the left, having an abnormally low [N ll] A6583/Ha ratio. This is because SBS 0335-052E 
has a low N/0 ab undance ratio, a characteristic of extremely low-metallicity BCDs (e.g. 
Izotov et al. II2006I ). and the ionizing parameter of its H ll region is high, resulting in a small 
abundance of singly ionized nitrogen. For comparison, the grey dots represent all galaxies 
from the SDSS DR7 with fiux errors smaller than 10% for each of the four emission lines 
H/9, [O III] A5007, Ha and [N ll] A6583. The emission line fiuxes for the SDSS galaxies 



shown by grey dots were measured using the technique developed by iTremonti et al. I (l2004j ) 



and are taken from the SDSS wetaiteH The.e gala.ie= are di.tribnted into two wing=, the 
left one thought to contain star-forming galaxies and the right one AGN hosts. The dashed 
lin e represents the eni p irical divisory line between star-forming galaxies and AGNs proposed 
by iKauffmann et al. I (|2003[ ) , wh ile the continuou s line represents the upper limit for pure 
star-forming galaxies taken from IStasihska et al. I (120061 ). We see that the LCGs lie in the 
low-metallicity part of what is usually considered as the re gion of star-formi n g gal axies. 
We note that the fi ve low-metallicity AGN candidates from llzotov fc Thuan I ( l2008l ) and 
Izotov et al. I (120101 ). represented by large filled circles, lie in the same region as the LCGs. 
The BPT diagram cannot distinguish between star-forming galaxies and AGN hosts in the 
low-metallicity regime. 



^http://www. sdss.org/DR7/products/value_added/index. html. 



- 9 - 



Element abundances 



We derived element abundances in all 803 LCGs from the narrow emission-line fluxes, 
using the so-called direct method. This method is based on the determination of the elec- 
tron temperature within the [O ill] zone from the [O lll]A4363/(A4959+A5007) line ratio. 
The fluxes in all spectra were measured using Gaussian fitting with the IRA F§| SPLOT 



routin e. They were corrected for both extinction, using the reddening curve of IWhitford 
(ll958l ). and underlying hy drogen stellar absorp tion, derived simultaneously by an itera- 



tive procedure described by llzotov et al. I ( 1l994j ) and using the observed decrements of the 
narrow hydrogen Balmer lines. The extinction coefficient C(H/3) and equivalent width of 
hydrogen absorption lines EW(abs) are derived in such a way to obtain the closest agree- 
ment between the extinction-corrected and theoretical recombination hydrogen emission-line 
fiuxes normalized to the H/3 fiux. It is assumed that EW(abs) is the same for all hydrogen 
lines. This assumption is justified by the evolutionary stellar population synthesis models of 



Gonzalez Delgado et al. I (|2005[ ) 



The physical conditions, and the ionic and total heavy element abundances in the H 



II regions were derived following llzotov et al. I ( l2006l ). In particular for O^"*", Ne^"*", and 
Ar^+ abundances, we adopt the temperature Te(0 ill) directly derived from the [O ill] 
A4363/(A4959 + A5007) emission-line ratio. Th e electron temperat ures Te(0 ll) and Te(S 



III) were derived from the empirical relations by llzotov et al. I ( l2006l ). We used Te(0 ll) for 
the calculation of 0+, N^, S^, and Fe^"*" abundances and Te(S ill) for the calculation of S^"*", 
CP+, and Ar^"*" abundances. The electron number densities A'"e(S ll) were obtained from the 
[S ll] A6717/A6731 emission-line ratios, respectively. The low-density limit holds for the H ll 
regions that exhibit the emission lines considered here. The element abundanc es then do not 



depen d sensitively on iVg. We use the ionization correction factors (ICFs) from llzotov et al. 



( I2OO6I ) to correct for unseen stages of ionization and to derive the total O, N, Ne, S, CI, Ar, 
and Fe abundances. 

The dependence of the oxygen abundance 12 +log 0/H of the LCGs on their redshift 
z is shown in Fig. There is no clear variation of 12-|-log 0/H with z. The median value 
of 12 -|- log 0/H for the whole sa mple is ~ 8.11 or 0.26 solar (Table [1]). This median value 



is similar t o the one obtained bv lAmorfn et al. I (|2010[ ) for the subsample of the green pea 



galaxies of ICardamone et al. I ( 120091 ) they used. However, our values of 12 + l og /H are 



significantly lower - by ~ 0.6 dex - than those obtained by ICardamone et al. I (|2009[ ). This 



^IRAF is the Image Reduction and Analysis Facility distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities for Research in Astronomy (AURA) under 
cooperative agreement with the National Science Foundation (NSF). 
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difference is due to t he different techniques us ed to derive abundances. We have used the 



direct method while ICardamone et al. I ( 120091 ) have used an empirical cahbration method 



based on the [N ll] A6583 /Ha flux ratio, despite the fact that nearly all their spectra contain 
a strong enough [O ill] A4363 emission line, which can be used to determine the electron 
temperature of the [O ill] zone. Evidently, there is an offset between oxygen abundances 



Cardamone et al. 


(20091). 


and by 


Cardamone et al. 



Fig. [6^ also shows that LCGs with lower equivalent widths EW(H/3) (open symbols) 
have systematically higher 12 + log 0/H than LCGs with higher equivalent widths (filled 
symbols). While the median 12 + log 0/H is 8.04 for subsample 1, it is 8.16 for subsample 
4 (Tabled]). This effect is more pronounced in Fig. where we show the dependence of 12 
+ log 0/H on EW(H/3). The LCGs with high EW(H/3) > 200A are more metal-poor than 
LCGs with lower EW(H/3) by ~ 0.2 dex. Such variations can be explained by the fact that 
LCGs with a lower EW(H/3) are in general more massive and thus more metal-rich (see the 
discussion of the mass-metallicity relation below). 

In Fig. d we show the dependence of the N/0, Ne/0, S/0, Cl/0, Ar/0 and Fe/0 
abundance ratios on the oxygen abundance 12 +log 0/H for the LCG sample (filled cir- 
cles). For comparison, we also show a sample of BCDs that has been collected primarily to 
study helium abundances in low-metallicity dwarf gal axies. This samp l e is th e same as in 
Izotov et al. I ( 119941 ). with the addition of galaxies from llzotov fc Thuan I ( 120041 ) and includes 
mainly nearby low-metallicity and low-luminosity galaxies. It is seen from Fig. [7] that the 
distribution of heavy element abundance ratios in LCGs is very similar to those in BCD 
galaxies, despite the fact that LCGs have higher oxygen abundances (12 + log 0/H > 7.6), 
while the BCD galaxies span a larger range of metallicities, with 12 + log 0/H as low as 
7.12. This implies that the chemical enrichment history has been very similar in BCDs and 
LCGs. 

The N/0 vs 0/H diagram (Fig. [7^) shows a significant spread of N/0 values, in the 
range -0.9 - -1.7 at 12 + log 0/H > 7.6, for both the LCG and BCD samples. However, 
there is no appreciable increase of N/0 with increasing 12 + log 0/H up to ~ 8.4, sug- 
gesting t hat nitrogen is r nainly of primary origin in LCGs. This is at variance with the 
results of lAmorm et al. I (120101 ) who find a clear increase of N/0 with 12 + log 0/H over 
the whole 7.5 - 8.4 range, suggesti ng that secondary production of nitrogen is important. 
Additi onallv. lAmorfn et al. I ( 120101 ) derived a high log N/0, larger than the solar ratio of 



-0.86 ( Asplund et al. 20091 ). in a significant number of green pea galaxies with 12 + log 



0/H in the range ~ 8.1 - 8.4. The reason for the differences between our results and those 
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of lAmorm et al. I ( 120101 ) is not clear. We note that lAmorm et al. I ( 120101 ) used the pipehne 
measurements of hne fluxes available in the SDSS archive, while we measured them directly 
using the SDSS spectra, taking into account the possible contamination of [N ll] A6583 by 
Ha A6563. 

It is worth noting the increase of the Ne/0 abundance ratio and the decrease of the 
Fe/0 ratio with increasing oxygen abundanc e. These trends are caused by the depletion of 



O and Fe onto dust grains, as discussed by llzotov et al. 
these trends. 



(120061 ). The LOG data confirm 



4. Luminosity-metallicity relation 



The relation between the oxygen abundance 12 + log 0/H and the extinction-corrected 
absolute magnitude Mg in the SDSS g band for LCGs is shown in Fig. |Si The symbols for 
LCGs are the same as in Fig. [2J The g apparent magnitudes are taken from the SDSS data 
base. The correction for extinction is done using the extinction coefficient C(H/3) obtained 
from the hydrogen Balmer decrement in the SDSS optical spectra. Distances of LCGs are 
derived from redshifts and adopting a Hubble constant Hn = 75 km s~^ Mpc~ ^. For compari- 
son, we show by dots the emission-line galaxy sample of iGuseva et al. I ( 120091 ). We also show 
by large filled circles the three most metal- deficient BCDs known in the local universe, I Zw 
18, SBS 0335-052W and SBS 0335-052E jGuseva et al. Iboooh . The intermedia t e-red shift 
{z < 1) extremely low-metallicity emission-line galaxies studied by lKakazu et al. I (j2007l ) are 
shown by l arge filled squares, w hile the luminous metal-poor star- forming galaxies at z ~ 0.7 



studied by iHoyos et al. I (120051 ) are shown by large filled triangles. Th e dotted-line rectangl e 



indicates the location of the Lyman-break galaxies (LBG) at z ~ 3 of iPettini et al. I ( l200ll ). 



The solid line is the luminosity-metallicity relation derived by lGuseva et al. 1 J2009h for their 
emission-line galaxy sample (solid line in their Fig. 9). 

LCGs are clearly offset by ~ 2 mag to brighter magnitudes as compared to the bulk of 
the local emission-line galaxies, implying strong ongoing star formation. There exists a sim- 
ilar offset for the three most-metal deficient BCDs I Zw 18, SBS 0335-052W and SBS 0335- 



052E, the intermediate-redshift extremely metal-poor emission-line galaxies of iKakazu et al. 



(120071 ). the l uminous metal - poor e mission-hne galaxies at z ^ 0.7 of iHoyos et al. 



(120051 ) and 



the LBGs of iPettini et al. I (120011 ). It appears that all these types of galaxies are undergoing 
strong bursts of star formation and that they define a common luminosity-metallicity rela- 
tion, the one for actively star-forming galaxies, going from the most-metal deficient BCDs 
at the faint end to LBGs at the bright end. The best linear likelihood fit to the strongly 
star-forming galaxies, excluding LBGs for which the oxygen abundances are derived by using 
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not the direct but the strong-hne method, is 

12 + logO/H= (-0.16 ±0.01) X + (4.86 ±0.19). (1) 

It is shown by a dashed hne in Fig. |8l 

Our LCGs, despite their lower redshift range (z ^ 0.02 - 0.6 3), are very similar in lumi- 



nosity to the luminous galaxies at z ~ 0.7 of iHoyos et al. I ( l2005l ). but they are, in the mean 



some 3 mag fainter than LBGs. However, the most luminous objects in our LCG sample 
do have oxygen abundances and luminosities that are similar to those of LBGs at 2; ~ 3, 
implying that these two types of objects are likely of the same nature. In particular, most 
of the LCGs inside the LEG dotted rectangle (Fig. |H]) have disturbed morphologies and rel- 
atively low equivalent widths EW(H/3). This implies some kind of gravitational interactions 
and/or merging processes in the galaxies with the significant fraction of a non- ionizing stellar 
population that contributes non-negligibly to the observed stellar continuum. The similarity 
in properties between LCGs and LBGs suggests that, since LCGs are relatively nearby, the 
bulk of them being in the redshift range - 0.3, they can be studied in detail to shed light 
on the physical properties of distant LBGs with ^ ~ 3. 



5. Stellar masses 

The stellar mass of a galaxy is one of its most important global characteristics. It 
can be derived by modeling its spectral energy distribution (SED) which depends on the 
adopted star formation history. However, in the case of LCGs, the situation is complicated 
by the presence of s trong ionized gas e missio n, which must be subtracted before stellar 



mass determination. ICardamone et al. I ( 120091 ) did consider the effect of the strong nebular 



emission lines. However, they did not take into account the effect of the gaseous continuum 
emission which is significant in spectra of LCGs with a high EW(H/3). Neglecting the 
correction for gaseous continuum emission would result in a significant overestimate of the 
galaxy stellar mass for two reasons: 1) gaseous continuum emission increases the luminosity 
of the galaxy; and 2) the SED of gaseous continuum emission is flatter than that of stars, 
making the SED redder than expected for pure stellar emission. Consequently, the fraction 
of the light from the red old stellar population is artificially increased. Therefore, to derive 
the correct stellar mass of the galaxy, we have used the method described below that takes 
into account the contribution of gaseous continuum emission. 
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5.1. Method 

The method consists of fitting; a series of model S EDs to the observed one and finding the 



best fit. It is described in lGuseva et al. I (120061 120071 ) and consists of the following. The fit is 



performed for each SDSS spectrum over the whole observed spectral range of AA3900-9200A, 
which includes the Balmer jump region (A3646A) for high-redshift LCGs with z > 0.1, and 
the Paschen jump region (A8207A) for LGCs with z < 0.12. As each SED is the sum of both 
stellar and ionized gas emission, its shape depends on the relative contribution of these two 
components. In LCGs with high EW(H/3) > lOOA (subsamples 1 and 2), the contribution of 
the ionized gas emission can be very large. However, the equivalent widths of the hydrogen 
emission lines never reach the theoretical values for pure gaseous emission, which for the 
H/3 emission line is ~ 900 - lOOOA, slightly depending on the electron temperature of the 
ionized gas. The objects in our LCG sample span a lower range of EW(H/3)s, between 50 and 
~500A (Fig. [6)d), due to the contribution of stellar continuum emission. The contribution of 
gaseous emission relative to stellar emission can be parameterized by the equivalent width 
of the H/3 emission line EW(H/3). Given a temperature Te(H"^), the ratio of the gaseous 
emission to the total emission is equal to the ratio of the observed EW(H/3) to the equivalent 
width of H/9 expected for pure gaseous emission. The shape of the spectrum depends also 
on reddening. The extinction coefficient for the ionized gas C(H/5)gas can be obtained from 
the observed hydrogen Balmer decrement. However, there is no direct way to derive the 
extinction coefficient C(H/3)stars for the stellar emission, which can be in principle different 
from the ionized gas extinction coefficient. Finally, the SED depends on the star formation 
history of the galaxy. 

We have carried out a series of Monte Carlo simulations to reproduce the SED of each 
galaxy in our sample. To calculate the contribution of the stellar e mission to the SEDs, 



we have adopted the grid of the Padua stellar evolution models by iGirardi et al. I (120001 ) 
with heavy element mass fractions Z = 0.001, 0.004, and 0.008. To reproduce the SED of 
the stellar c omponent with any star formatio n history, we have calculated with the package 



PEGASE.2 JPioc fc Rocca-Volmerang^ll997h a grid of instantaneous burst SEDs for a stellar 



mass of 1 Mq in a wide range of ages, from 0.5 Myr to 15 Gyr. We have adopted a stellar 
initial mass function with a Salpeter slope, an upper mass limit of 100 Mq, and a lower 
mass limit of 0.1 Mq. Then the SED with any star formation history can be obtained 
by integrating the instantaneous burst SEDs over time with a specified time-varying star 
formation rate. 

We have approximated the star formation history in each LCG by a recent short burst 
with age t (young) < 10 Myr, which accounts for the young stellar population, and a prior 
continuous star formation responsible for the older stars, with age starting at t2 = t(old) 
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and ending at ti, with t2 > ti and varying between 10 Myr and 15 Gyr. The contribution 
of each stellar population to the SED is parameterized by the ratio of the masses of the old 
to young stellar populations, b = M (old) /M (young), which we vary between 0.01 and 1000. 

The contribution of the gaseous emission was scaled to the stellar emission using the ratio 
of the observed equivalent width of the H/3 emission line to the equivalent width of H/3 ex- 
pected for pure gaseous emission. Then the total modeled monochromatic (gaseous+stellar) 
continuum flux near the H/9 emission line for a mass of 1 Mq is scaled to fit the monochro- 
matic extinction-corrected luminosity of the galaxy at the same wavelength. The scaling fac- 
tor is equal to the total stellar mass M* in solar units. In our fitting model M*= M(young) 
-|- M(old), where M(young) and M(old) are respectively the masses of the young and old 
stellar populations in solar units. These masses are derived using and b. 



The SED of the gaseous continuum is taken from lAUer I (Il984j ) and includes hydrogen 
and helium free-bound, free-free, and two-photon emission. In our models, it is always 
calculated with the electron temperature Te(H"'") of the H"*" zone and with the chemical 
composition derived from the H ll region spectrum. The observed emission lines corrected 
for reddening and scaled using the flux of the H/3 emission line were added to the calculated 
gaseous continuum. The flux ratio of the gaseous continuum to the total continuum depends 
on the adopted electron temperature Te(H"'") in the H+ zone, since EW(H/3) for pure gaseous 
emission decreases with increasing Te(II+). Given that Te(H"'") is not necessarily equal to 
Te(0 III), we chose to vary it in the range (0.7-1.3) xTe(0 ill). Strong emission lines in 
LCGs (see Fig. [T]) are measured with good accuracy, so the equivalent width of the H/3 
emission line and the extinction coefficient for the ionized gas are accurate to 5% and 20%, 
respectively. Thus, we vary EW(H/3) between 0.95 and 1.05 times its nominal value. As for 
the extinction, we assume that the extinction coefficient C(H/3)stars for the stellar light is the 
same as C(H/3)gas for the ionized gas, and we vary both in the range (0.8-1.2) xC(H/3), where 
C(H/3) is the extinction coefficient derived from the observed hydrogen Balmer decrement. 
For each LCG, we run 10^ Monte Carlo models varying t(young), ti, t(old), b, and Te(II+) 
randomly in a large range, and EW(H/9) and C(H/9)gas in a relatively smaller range because 
the latter quantities are more directly constrained by observations. The best modeled SED 
is found from minimization of the deviation between the modeled and the observed 
continuum in five wavelength ranges which are free of the emission lines. 

In Fig. [9^ we show the SED fitting of the spectrum of the galaxy J085 1+5840 with 
EW(H/3) = 303A. The total modeled spectrum is labeled "total" while the stellar and gaseous 
spectra are labeled "stars" and "gas" , respectively. The model SED fits very well the observed 
spectrum of the galaxy. It is seen that the contribution of gaseous emission is high because 
of the high EW(H/3). The fraction of the gaseous emission as a function of wavelength in the 
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spectrum of the same galaxy is shown in Fig. [Qjo. The two steep jumps seen at ~ A3660A 
and ~ A8200A are respectively due to the Balmer and Paschen discontinuities of the ionized 
gas emission. Between these jumps, the fraction of gaseous emission increases monotonously 
from ~ 20% to ~ 50% making the spectrum flatter than a pure stellar emission spectrum. 
Neglecting the contribution of this gaseous emission would have resulted in an artificial 
increase of the galaxy mass. For the stellar mass determination, we use the SED labeled 
"stars" in Fig. [9^, which is considerably below the SED labeled "total" . 

In Fig. [TUl we show the dependence of (a) the absolute g magnitude Mg and (b) of 
the total stellar mass on the fraction of gaseous continuum near the H/5 emission line. 
The solid lines are the best linear likelihood fits to the data. It is seen that galaxies with 
high EW(H/9) (filled symbols) are on average fainter and less massive than galaxies with low 
EW(H/9) (open symbols). 



5.2. Masses of young and old stellar populations 



From the SED fits of the spectra of all the 803 LCGs in our sample, we find the median 
total stellar mass to be M* = 1x10^ Mp) (Table [H). This value is a factor of ~ 3 lower than 
the median mass 3x10^ Mq found by lCardamone et al. I ( 120091 ) for their sample of green pea 
galaxies. As noted before, these star-forming green pea galaxies with their round morphology 
and high EW(H/3) are similar to the LCGs in our subsample 1 which is characterized by 



an even smaller median m ass of 7x10^ Mp^ 



median mass between the ICardamone et al. I f 



Tabl e [T]). We believe the factor of ~4 in 
2009!) g r een p ea sample and our subsample 



1 is primarily due to the fact that ICardamone et al. I (120091 ) did not subtract the gaseous 
continuum emission from the total galaxy emission in deriving th eir stellar masses. A more 
detailed comparison between stellar masses derived by us and by ICardamone et al. I (120091 ) 
is given in Section 5.3. 

As for the young stellar population in LCGs, its median mass for the whole sample 
is M(young) = 3x10^ M©, i.e. a factor of ~ 30 below the median mass of the old stellar 
population. There is a small increase of the median mass of the young stellar population by 
a factor of ~ 1.4 from subsamples 1 and 2 to subsamples 3 and 4 (Tabled]). For comparison, 
the median total stellar mass of LCGs in subsample 1 is ~ 2 times smaller than that in 
subsample 4. 

Fig. [TTh shows the dependence of the total stellar mass M^, of LCGs on the mass 
M(young) of the young stellar population. We note that the total stellar mass is in the 
range 10^ - 5x10^*^ Mq, while the mass of the young stellar population is in the range 2x10^ 
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- 10^ Mq. Evidently, there is an increase of the mass of the young stellar population with 
increasing total stellar mass of the galaxy. This suggests that more massive galaxies contain 
more regions of ongoing star formation. Fig. [TTb shows the dependence of the total stellar 
mass on the ratio b of the mass of the old stellar population to that of the young stellar 
population. It is seen that M (old) /M (young) increases with increasing M*, reaching the 
high values of several hundreds, predominantly in galaxies with low EW(H/3) (open symbols). 
On the other hand, there is an appreciable number of LCGs with high EW(H/3) and low 
M(old)/M (young) (filled symbols), implying that less massive dwarf galaxies {10'^ Mq < 
< 10^ Mq) are also less evolved, both in terms of their stellar populations and metallicities. 

In the fitting analysis and in the determination of global characteristics of LCGs such as 
luminosities and masses, we have made no correction for aperture effects. Yet, these effects 
can be important, especially for galaxies at low redshifts, which may have angular sizes 
larger than the 3" diameter of the fibers with which the SDSS spectra were obtained. We 
estimate these aperture corrections below. Assuming that the spectral energy distribution 
in the galaxy inside and outside the round aperture are the same, the aperture correction 
factor can be written as 

fiapcorr) = iQO-^CKs'VKtot))^ ^2) 

where r(3") and r(tot) are respectively the SDSS r apparent magnitude inside the 3" round 
aperture and the modeled integrated SDSS r apparent magnitude of the galaxy as derived 
from the SDSS images. Both r(3") and r(tot) are taken from the SDSS DR7 database. Using 
Eq. [21 we found the median flux correction factor f{apcorr) to vary from ~ 1.4 for subsample 
1 to ~ 1.8 for subsample 4. However, for comparison purposes, we have decided not to apply 
these aperture corrections to the LCG sample, as they are not taken into account in other 
studies. Furthermore, the assumption that the spectral energy distributions of the galaxy 
inside and outside the slit are the same may not be valid. Thus, we should keep in mind 
that the true global characteristics for LCGs, such as luminosity, mass or star formation 
rate, may be larger than those derived here from the galaxy spectra, especially for nearby 
galaxies. 



5.3. Comparison of stellar masses and o xygen abundances o f gree n pea 



galaxies derived by us and by ICardamone et aL I (120091 ) 



As discussed before, there are 66 green pea galaxies from ICardamone et al. I (|2009[ ) 



that are included in the LCG sample. We c ompare here their stellar masses and oxygen 
abundances derived by us to those obtained by ICardamone et al. I (120091 ). We show in Table 
|2] the SDSS ID number of each galaxy, its equatorial coordinates, its stellar masses and 
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oxygen abundances derived by lCardamone et al. I (120091 ) and in this paper. An examination 
of the Table shows that the stellar masses o f green pea galaxie s deriy ed in this paper are 



systematically lower than those obtained by ICardamone et al. I (120091 ). This effect is also 



clearly seen in Fig. [T^ . As already discussed, this systematic offset is due to the fact that, in 
fitting the SED, we subtract the contribution of gaseous continuum emission, which results 
in lower stellar masses. This effect is largest for lower-mass galaxies with log M^jM^ < 8.5. 
These are mostly galaxies with high EW(H/3) in subsample 1. 

Fig. W2o show s the c omparison of oxygen abundances for the same galaxies. Again, the 



Cardamone et al. I (120091 ) abundances are systematically larger (~ 0.6 dex in 12-|-logO/H) 
than ours. This systematic offset is a consequence of the different abundanc e deter mination 
methods used, the more uncertain strong-line method by lCardamone et al. I (|2009[ ). and the 



more reliable direct method by us. We note that nearly all abundance determinations in 
high-redshift galaxies and in the majority of nearby SDSS star-forming galaxies are based on 
strong-line methods. The more uncertain oxygen abundances that result from these methods 
should not be used in galaxy evolution studies without a careful analysis of their errors. 



6. Younger starbursts are hotter and denser than older ones 

We investigate here how the properties of the luminous H ll regions in LCGs depend 
on the stellar mass of the galaxy. In particular, we examine the dependence of the electron 
temperature re(0 ill) and the electron number density A'"e(S ll) on the total stellar mass 
of the galaxy (respectively Figs. [T3k and [T3b). We find that both the electron temperature 
and the electron number density are systematically higher in galaxies with higher EW(H/9). 
These galaxies also have systematically lower masses. The median electron temperatures 
and electron number densities of the galaxies in the different LCG subsamples vary from 
13 lOOK and 180 cm~^ in subsample 1, to 11 700K and 90 cm~^ in subsample 4. 

Young stellar populations in galaxies with higher EW(H/3) have younger ages t(young) 
(Tabled]). Thus, we conclude that the younger starbursts in our LCGs are systematically 
hotter and denser. 



7. Mass-metallicity relation 

In Fig. [T3^,b we show the mass-metallicity relation for LCGs. It is seen that galaxies 
with lower mass have systematically lower oxygen abundances. These lower-mass galaxies 
are primarily from subsamples 1 and 2 with high EW(H/3), while LCGs with lower EW(H/3) 
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are higher-mass objects with higher oxygen abundance. Although there is a large scatter , 
the data can be fitted by a linear relation, shown in Fig. [T^ by a solid line 

12 + logO/H= (7.475 ± 0.426) + (0.069 ± 0.047) x logM,, (3) 

where is in solar units. We also show by a dashed line in Figure [Tik the polynomial 



fit obt ained by I Amor in et al. I ( l2010l ) to the green pea galaxy sarnple of ICardamone et al. 



( I2OO9I ). The difference between our fit and that of lAmorm et al. I (120101 ) is due to the fact 
that we take into account the gaseous continuum emission and thus derived lower stellar 
masses produc i ng ari offset by a factor of ~ 4 in the galaxy mass between our data and 



Amorm et al. I (|2010[ ) fit. Furthermore, because of a small stellar mass ra nge and a large 



disper sion in oxygen abundances, there is no strong evidence from Fig. 3 of lAmorm et al. 



( I2OIOI ) t hat the data for green pea galaxies is better fit by a polynomial than a line. A linear 



fit to the lAmorm et al. I (|2010[ ) data would result in a better agreement with the fit to LCGs 



derived in this paper. 



We show by a dash-dotted line the fit of lAmorm et al. I (120101 ) to a sample of star- 
forming galaxies assembled from the SDSS. This relation is considerably steeper than the 
one obtained for LCGs. The difference is probably not due to the determination of the stellar 
masses of the SDSS star-forming galaxies. The equivalent width of their H/3 emission line is 
generally small, so that ionized gas emission is small. However, the oxygen abundances of 
these galaxies are derived using strong-line empirical methods which produce a systematic 
offset as compared to the oxygen abundances derived by the direct method. That is probably 
why the relations for SDSS star-forming galaxies and LCGs differ. 

The relation shown by the solid line in Fig. [T^ is also at odds with the finding of 
Cardamone et al. I ( 120091 ) for their green pea sample. Those authors did not find a clear 
evidence for the dependence of the oxygen abundance on stellar mass. This is likely be- 
cause of the way they derive the metallicities and masses of their galaxies. As discussed 
before, oxygen abundances were derived by using the more uncertain strong-line technique, 
and stellar masses were determined without subtracting the gaseous continuum emission 
from the observed SED. In that way, both oxygen abundances and galaxy stellar masses 
are overestimated. Ano ther reason for why a mass-metallicity relation was not seen by 
Cardamone et al. I ( 120091 ) is probably because the 10^ - 10^°'^ Mq mass range of their sample 
is smaller than the 10'' - lO^'''^ Mq mass range of our LCG sample. 

In contrast to Fig. Ek, Fig. \\M ) shows that the N/0 abundance ratio for LCGs with 
EW(H/3) > lOOA (filled symbols) is independent on stellar mass, implying a primary origin 
for nitrogen in these galaxies. As for LCGs with EW(H/3) < lOOA (open symbols), there 
is a dependence of N/0 ratio on the total galaxy mass, implying some secondary nitrogen 
production in the highest-mass galaxies. 
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8. Star formation rate 



One of the most important characteristics governing a galaxy's evolution is its star 
formation rate. We derive the star formation rate SFR(Ha), using the extinction-corrected 
luminosity L(Ha) of the Ha emission line and the star formation rate SFR(FUV ), using the 



extinc tion-corrected GALEX FUV luminosity L(FUV), and relations given by iKennicutt 
(|l998h 



SFR(Ha) 
SFR(FUV) 



7.9 X lO-^'L(Ha), 
1.4 X 10"^^L(FUV). 



(4) 
(5) 



In equations H] and 0, SFR is in units of Mq yr ^, L(Ha) in erg s ^ and LfFUV) in erg 
s~^ H z~^. To correct L(FUV) for extinction we adopt the reddening law of ICardelli et al. 
Jl989h with R = A{V)/E{B -V) = 3.1 and the extinction coefficient C(R/3) derived from 



the optical SDSS spectrum. Then the extinction correction factor for L(FUV) is equal to 
2.512^(™^) with A(FUV) = 8.24x^(5 - V) = 5.58xC(H/3). 

The GALEX FUV fluxes are available for 675 galaxies out of the 803 LCGs in our sample. 
In Table [H we give the median SFR(Ha) and SFR(FUV) for the different subsamples and 
for the total sample. It can be seen that the SFR(Ha) and SFR(FUV) are in good agreement 
and that their median values do not change appreciably from one subsample to the next, 
being ~ 4 Mq yr~^. Because of the good agreement, we will consider below only SFRs 
derived from the Ha luminosity. 

In Fig. [T5k we show the dependence of SFR(Ha) on the total stellar mass M^,. The 
star formation rate for LCGs is in the range 0.7 - 60 Mq. The median values in different 
subsamples are in the narrow range 3.4 - 5.0 Mq yr~^ (Tabled]). Thus, there is no significant 
difference in the average SFR for galaxies with different EW(H/3). However, we note the 
increase of SFR(Ha) with increasing total stellar mass in galaxies from subsamples 3 and 4 
(galaxies with EW(H/3) < lOOA shown in Fig. [T5b by open symbols), while there is a less 
strong dependence for galaxies from subsamples 1 and 2 (galaxies with EW(H/3) > lOOA 
shown in Fig. [TBb by ffiled symbols). The ran ge of SFRs in LCGs is comparable to th at in 
the interr nediate-redshift star- forming galaxies (IHoyos et al. Il2005t iKakazu et al. 1120071 ) and 



in LBGs ( iPettini et al. Il200ll ). but is ~ 10-100 times higher than in typical nearby BCDs. 
For example, the SFR in the prototype BCD I Zw 18 is 0.1 Mq yr'^ jThuan hom . 



The dependence of the specific star formation rate SSFR(Ha) = SFR(Ha)/M* is shown 
in Fig. [T5b . The SSFR is of interest because it is related to the time needed to double the 
stellar mass of a galaxy. We show by dashed lines from left to right the loci of constant 
star formation corresponding to 1, 10 and 100 Mq yr~^, respectively. It is seen that there 
is an inverse dependence of the SSFR on the total stellar mass of LCGs. The lowest SSFRs 
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(~ 10"^ yi'^^) occur in the highest-mass galaxies with EW(H/3) < lOOA, while the highest 
SSFRs (up to ~ 10~^ yr~^) are found almost exclusively in galaxies with EWfH/j ) > lOOA. 
These values are very high, higher than those found by ICardamone et al. I (120091 ) for their 
green pea galaxies (their maximum SSFRs are ~ 10~^ yr~^), mainly because those authors 
overestimated the galaxy stellar masses. They are also considerably higher than those found 
in merging galaxies which lie in the 10~^° yr~^ - 10~^ yr~^ range. Our SSF R values are 



comparable to those found in high-redshift galaxies in the ranges z = 4-6 (jStark et al. 



20091 ) and 6-8 (jSchaerer fc de Barros Il2010l ). Thus, our LCGs resemble galaxies with the 
highest SSFR known. 

The quantity l/SSFR(Ha) is also a measure of the time needed to double the metallicity 
of a galaxy. Thus, if a constant star formation rate is assumed, then the ISM in our galaxies 
is enriched by heavy elements on a time scale of ~ 10^ - 10^ yr, considerably smaller than 
the Hubble time. It is likely then that star formation in LCGs proceeds in short bursts of 
duration 10^ - 10^ yr. Then, the SSFR averaged over time would be much lower. Bursts of 
star format ion separated by long quiescent periods of several Gyr (not unlike the situation in 
BCDs, see iThuan 1120081 ) are also required in the LCGs with the highest SFRs, so as not to 
exhaust the gas supply for star formation. Assuming continuous star formation, then with 
a hydrogen gas supply of 10^° Mq and a SFR of 50 Mq yr~^, the time for gas depletion is 
only ~ 2x10^ yr. Clearly a burst mode is required for the LCGs with the highest SFRs. 



9. Age of stellar populations and the burst mode of star formation in LCGs 

While the young stellar population in our galaxies is modeled by a short burst occurring 
at i:(young) ~ 3 - 4 Myr (Table [T]), the older stellar population is fitted by a model with 
continuous star formation starting at t2 = t(old) and ending at ti, with t2 > ti. The 
parameter t(old) sheds light on the evolutionary status of the galaxy (whether the galaxy 
is young or old), and helps to understand the character of its star formation (whether it 
proceeds in a bursting or continuous manner). 

Fig. [T6] shows how the age of the oldest stellar population )f:(old) depends on the total 
stellar mass M^,. The dashed lines are the loci of constant continuous star formation rate 
corresponding to, from left to right, SFR = 1, 10 and 100 Mq yr~^. 

It can be seen that galaxies with high EW(H/3) (> lOOA) denoted by filled symbols 
do not show any evident trend. On the other hand, there is an increase of t(old) with 
increasing for galaxies with EW(H/3) < lOOA. This implies that, among the galaxies 
with low EW(H/3), massive galaxies are in general older than lower mass galaxies. The 
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absence of such a dependence for galaxies with high EW(H/3) may in part be explained by 
the uncertainties in the SED modeling. Indeed, the fraction of the light produced by the 
old stellar population in spectra of galaxies with a high EW(H/3) is much lower than that 
in galaxies with a low EW(H/9), both because of the brighter young stellar population and 
of the higher contribution of gaseous emission. Furthermore, galaxies with high EW(H/3) 
have systematically higher redshifts compared to galaxies with low EW(H/3). Their larger 
distances make their continua weaker, which renders their SED fitting more uncertain. This 
is because of the following selection effect. We select only galaxies with [O ill] A4363 that 
can be measured with good accuracy for abundance determination, and this line is stronger 
in galaxies with high EW(H/3). 

The location of the galaxies in Fig. [12] relative to the dashed lines indicates continuous 
SFRs that are somewhat lower than the ongoing SFRs derived by equations H] and 0, and 
shown in Fig. [T5b . The largest offset is for low-mass galaxies with high EW(H/3) (filled 
symbols). The vast majority of these galaxies are located in the region corresponding to 
continuous SFRs < 1 M© yr~^, i.e. > 10 times lower than those derived from the ongoing 
star formation episode. This discrepancy is again an indication that star formation in the 
LCGs occurs in rare bursts of short duration, so that the SFR averaged over the lifetime 
of the galaxy is more than one order of magnitude lower than the one derived from the 
Ha luminosity. The bursting nature of star formation is less pronounced for the higher- 
mass galaxies with low EW(H/3) (open symbols): their continuous SFRs averaged over the 
galaxy's lifetime are only ~ 2 - 3 times lower than those derived from the ongoing episode 
of star formation. 



10. Conclusions 



We discuss here the properties of a sample of luminous compact star-forming galaxies 
(LCG) selected from the Data Release 7 (DR7) of the Sloan Digital Sky Survey (SDSS). We 
show that the so-called "green pea" galaxies discussed by Cardaraone et al. (j2009 ) are a 
subset of these LCGs. However, in contrast to Cardamone et al. ( 20091 ) who selected their 
green pea galaxies solely on the basis of broad-band SDSS colors and compact structure, 
our LCG galaxies were chosen also on the basis of their spectroscopic properties. Our main 
selection criteria are as follows: a) a high extinction-corrected luminosity L{B./3) > 3x10^° 
erg s~^ of the H/3 emission line; b) a high equivalent width EW(H/3) > 50A of the H/9 emission 
line; c) a strong [O ill] A4363 emission line allowing accurate abundance determination; d) 
a compact structure on SDSS images; and e) an absence of obvious AGN spectroscopic 
features. Using these criteria, we selected in total 803 galaxies, split into four subsamples 
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based on the value of EW(H/9) (greater or smaller than lOOA) and morphological appearance 
(round shape or some evidence of disturbed morphology). 

Analysis of the LCG sample has led to the following findings: 

1. The selected galaxies have redshifts between 0.02 and 0.63, a redshift range that 
is a factor > 2 highe r thari the redshift range z = 0.112 - 0.360 of the green pea sample 
of ICardamone et al. I ( 120091 ). We find that the properties of LCGs and green pea galaxies 
selected from the SDSS are in many respects similar to those of well-studied nearby luminous 
blue compact dwarf (BCD) galaxies, such as SBS 0335-052E and HS 0837+4717. Since the 
color of an emission-line galaxy on the gri SDSS composite image changes with its redshift, 
depending on which strong emission lines fall within a particular filter band, the LCGs 
display a variety of colors from intense blue for nearby galaxies with z < 0.07, through pink 
and green at redshifts z ~ 0.09 and ~ 0.15 - 0.30, respectively, to brown at redshifts z ~ 0.6. 
Therefore, the green pea galaxy sample is just a subsample with a narrow redshift range of 
a larger LCG sample. A significant fraction of our LCGs have low redshifts {z < 0.1), thus 
allowing the study of LC G pr operties in mor e deta il than in the previous investigations of 
Cardamone et al. I (120091 ) and lAmorm et al. I ( 120101 ) . 



2. Although our selection criteria do not set an upper limit on the luminosity L(H/3) of 
the H/3 emission line, we find that there is an upper bound to L(H/3) at ~ 3x10^^ erg s~^. 
This implies a self-regulating mechanism of star formation that does not allow the formation 
of star-forming regions brighter than that limit in a galaxy. 

3. In the [O ill] A5007/H/3 vs [N ll] A6583/Ha diagnostic diagram, LCGs occupy the 
region of star-forming galaxies with the highe s t excitatiori . On the other hand, the five low- 
metallicity AGNs discovered by llzotov et al. I ( 120071 . |2010| ) and llzotov fc Thuan I ( l2008l ) also 
lie in the same region. This implies that further detailed studies of LCGs may result in the 
discovery of a larger number of low-metallicity AGNs showing some Seyfert 1 features such 
as broad hydrogen emission lines. 

4. We find that the oxygen abundances 12 + log 0/H in L CGs are in the range 
7.6 - 8.4 with a median value of ~ 8.11, confir ming the results of 



concerning a subset of the green pea sample of Cardamone et al. ( 20091 ). This range of 



Amorm et al. I (120101 ) 



oxygen abundances is typical of nearby lower-lum inosity BCDs. We note, h owever, that the 



median oxygen abundance of ~ 8.7, derived by ICardamone et al. I ( l2009l ) for their green 



pea sample is overestimated, because those authors used empirical r elations for abundanc e 
determinations, instead of the direct method used in this paper and in lAmorm et al. I ( 120101 ). 
There is no dependence of 12 + logO/H on redshift. On the other hand, we find that the 
oxygen abundances in galaxies with higher equivalent width EW(H/3) are generally lower. 
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As for the N/0, Ne/0, S/0, Cl/0, Ar/0 and Fe/0 abundance ratios, they are similar to 
those found for lower-luminosity BCDs. Overall, we find no appreciable difference in element 
abundances between LCGs and BCDs, implying a similar chemical enrichment history. 

5. In the luminosity- metallicity diagram, LCGs are shifted by ~ 2 mag to brighter 
magnitudes as compared to the bulk of nearby emission-line galaxies. However, together with 
the most metal-deficient nearby BCDs I Zw 18, SB S 0335-052W, SBS 03 35-052E, extremely 



metal-poor intermediate-redshi ft galaxies 
me tal-poor galaxi es (z ~ 0.7 



z < 1, 



Kakazu et al. 1120071 ). luminous compact 
Hoyos et al. II2005I ) and Lyman-break galaxies (LBG) {z ~ 



3, iPettini et al. Il200l[ ). LCGs form a common luminosity-metallicity relation, that for the 



most actively star-forming galaxies. Some LCGs have oxygen abundances and luminosities 
similar to LBGs, despite much lower redshifts. This opens the possibility of studying LBGs 
through brighter and hence easier to observe LCGs. 

6. We develop a technique for determing galaxy masses by fitting the observed spectral 
energy distribution (SED) by a model SED that includes both emission from stellar popu- 
lations of different ages and the continuum emission from the ionized gas, the latter having 
not been considered by previous investigators. We show that such an omission results in 
an overestimate of the stellar mass. This effect is especially strong for galaxies with high 
EW(H/3) > lOOA. Thus, the median mass of the gala x ies iii our sample i s a factor of 



3 time s lower than that derived by ICardamone et al. 
J2010h . 



(120091 ) and used by lAmorm et al 



7. The total stellar mass of LCGs is in the range ~ lO''' - ~ 10^^ Mq, with a median 
value of ~ 10^ Mq. The median mass of ~ 3x10'' Mq of the young stellar population is a 
factor of ~ 30 lower. We find that the mass of the young stellar population increases with 
increasing total stellar mass. However, the fraction by mass of the young stellar population 
is higher in galaxies with a low total stellar mass. Furthermore, H ll regions in low-mass 
LCGs are systematically hotter and denser, than those in high-mass galaxies. 

8. We find a clear increase of the oxygen abundance with increasing total stellar mass 
of the galaxy in the mass-metallicity diagram. However, there is no correlation between the 
N/0 abundance ratio and the total stellar mass for the lowest-mass galaxies, those with 
EW(H/3) > lOOA, implying primary N production in those objects. There is however a 
total mass dependence of N/0 for the highest-mass galaxies, those with EW(H/3) < lOOA, 
implying some secondary of N in those galaxies. 

9. The star formation rate in LCGs varies in a large range ~ 0.7 - ~ 60 Mq yr^^ with a 
median value of ~ 4 Mq yr^^. The range and median value of the SFRs i n LCGs is lower by a 
factor of several, as compared to Lyman-break galaxies at redshift ~ 3 ( IPettini et al. Il200l[ ) 
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and higher- redshift galaxies (IStark et al. 1 120091 : ISchaerer &: de Barros Il2010[ ). However, the 



highest SFRs found in LCGs compare well with those in high-redshift galaxies. LCGs have 
specific star formation rates that are among the highest known, in the range ~ 10~^ - 10~^ 
yr~^, and comparable to those of high-redshift galaxies in the 2; ~ 3 - 8 redshift range. This 
similarity offers the possibility of studying extreme star formation in high-redshift galaxies 
by using LCGs as relatively nearby proxies. 

10. We find that the age of the oldest stellar population increases with increasing stellar 
mass in the LCGs with low EW(H/3) (< lOOA). However, no such dependence is found for 
galaxies with high EW(H/3), an absence which can be explained by larger uncertainties in 
fitting the spectra of these galaxies. Star formation in LCGs at the rate derived from Ha 
emission cannot be continuous, but rather must occur in strong short bursts separated by 
long quiescent phases, just as in BCDs. The bursting nature of star formation is more 
pronounced in LCGs with higher EW(H/3) and lower masses. 
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Table 1. Average characteristics of LCGs 
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Table 2. Characteristics of ICardamone et al. I ( 120091 ) Green Pea galaxies included in the 

LCG sample 
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Table 2 — Continued 
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Fig. 1. — SDSS spectrum of a luminous compact emission- line galaxy in our sample obtained 
with the 2.5m Apache Point Observatory Sloan telescope in a Ih exposure. The ordinate 
scale is for the lower spectrum. The upper spectrum is the lower spectrum expanded by a 
factor of 100 along the ordinate. Weak emission lines with fluxes of ~ 1% that of H/3 are 
clearly detected in this spectrum. 
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Fig. 2. — Extinction-corrected H/3 luminosity I/(H/3) vs. redshift z for the LCG sample. 
LCGs with EW(H/9) > lOOA are shown by filled symbols. Those with a round shape and 
those with some signs of disturbed morphology are shown respectively by small filled cir- 
cles and small filled triangles. LCGs with 50A < EW(H/3) < lOOA are shown by open 
symbols. Those with a round shape and those with some signs of disturbed morphology 
are shown respectively by small open circles and small open triangles. Well-studied local 
LCGs, also known as Blue Compact Dwarf galaxies, such as Mrk 930, SBS 0335-052E, HS 
0837 +4717 and II Zw 40 are sho wn by large op e n circ les. Low-metallicity AGN candidates 
from llzotov fc Thuan I ( l2008l ) and llzotov et al. I (120101 ) are shown by large filled circles. 
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Fig. 3. — 50"x50" SDSS images of LCGs in order of increasing redshift from left to right. 
LCGs from subsample 1 with a round shape are shown in the upper panel, while LCGs from 
subsample 2 with some evidence of elongated structure and/or disturbed morphology are 
shown in the lower panel. 
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Fig. 4. — ((yf — r) - (r — color-color diagram for LCGs. The symbols for the various 
ty pes of objects are the sam e as in Fig. |2l Galaxies common with the green pea sample 



fcoosh 



of ICardamone et al. I (l2009[ l are encircled. Low-metallicity AGNs from llzotov fc Thuan 



are shown by larg e filled circles. The two solid lines are the dividing lines used by 



Cardamone et al. I ( l2009l ) for selection of green pea galaxies. For comparison, a repre sent a- 
tive sample of SDSS "normal" galaxies (black dots) and QSOs ( ISchneider et al. Il2010l ) (grey 
dots) is also shown. 
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Fig. 5. — The Baldwin-Phillips- Terlevich (BPT) diagram ( iBaldwin et al. Ill98ll ) for narrow 
emission lines. LCGs are shown by small open and filled symbols, their meaning being the 
same as in Fig. El Also, plotted are the 100 000 er aission-line galaxies frora SPS S DR7 (cloud 
of grey dots), the five low-met allicity AGNs from llzotov &: Thuan I (120081 ) and llzotov et al. 
J2010h (filled circles), and the two well-st u died B CDs SBS 0335-052E an d HS 0837+4717 



The d ashed line from iKauffmann et al. I (|2003[ ) and the solid line from IStasihska et al 



(120061 ) separate star-forming galaxies from active galactic nuclei. 
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Fig. 6. — Oxygen abundance 12 + log 0/H (a) vs. redshift z and (b) vs. the equivalent 
width EW(H/3) of the H/3 emission line. The LCGs are shown by symbols with the same 
meaning as in Fig. |2j 
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Fig. 7.— Abundance ratios (a) log N/0, (b) log Ne/0, (c) log S/0, (d) log Cl/0, (e) log 
Ar/0 and (f) log Fe/0 vs. oxygen abundance 12 + log 0/H for emission-line galaxies. Large 
open circles denote galaxies from the HeBCD sample (jlzotov et al. Ill994j ; llzotov fc Thuan 
20041 ). and small filled circles represent LCGs (this paper). 
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Fig. 8. — Luminosity-metallicity relation for emission-line galaxies. LCGs are shown by 
the same symbols as in Fig. [2] (filled a nd open circles and triangles). The emission- line 



galaxies studied by iGuseva et al. I (|2009[ ) are sh own by dots. The lu minosity-metallicity re- 
lation for these galaxies is fitted by a solid line (IGuseva et al. II2009I ). The most metal-poor 
local BCDs - I Zw 18, S BS 0335-052W and SBS 0335-052E- are shown by large filled cir- 



cles ( Guseva et al. IboOOl ). The extremely metal -poor emission- line galaxies at intermediate 
redshifts z < 1 studied by iKakazu et al. I (120071) an d the luminous metal-poor star-forming 



galaxies at z ~ 0.7 studied by iHoyos et al 



( 120051 ) are shown by large filled squares and 



large filled triangl es, respectively . The region occupied by Lyman-break galaxies (LBGs) at 
z ^ 3 studied by IPettini et al. I (1200 ll ) is indicated by a dotted rectangle. The linear best 
likelihood fit to the strongly star-forming galaxies, excluding LBGs, is shown by a dashed 
line. 
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Fig. 9. — (a) Best-fit model SED (solid black line labeled "total") superposed on the redshift- 
and extinction-corrected spectrum of the LCG SDSS J085H-5840 (grey line). The separate 
contributions from the stellar and ionized gas components are shown by black solid lines and 
labeled "stars" and "gas", respectively, (b) Gaseous emission fraction vs. wavelength for 
the modeled spectrum of LCG SDSS J0851+5840 with EW(H/3) = 303A. The two jumps 
seen at ~ A3660A and ~ A8200A are due respectively to the hydrogen Balmer and Paschen 
discontinuities in the ionized gas emission. 
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Fig. 10. — (a) The absolute magnitude Mg vs. the fraction of gaseous continuum near the 
H/3 emission line, (b) The total stellar mass vs. the fraction of gaseous continuum near the 
H/3 emission line. In both panels, the straight line represents the linear best likelihood fit. 
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— (a) The total stellar mass vs. the young stellar population mass M(young). 
total stellar mass vs. the young-to-old stellar population mass ratio. The masses of 
galaxies are dominated by an old stellar population, while the masses of low-mass 
are dominated by a young stellar population. 
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Fig. 12. — Com parison of: (a) stel l ar ma sses and (b) oxygen abundances derived for green 
pea galaxies by ICardamone et al. I ( l2009l ) and in this paper. Solid hues in both panels are 
lines of equal values. 
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Fig. 13. — The total stellar mass vs. (a) the electron temperature Te(0 ill) and (b) 
the electron number density iVe(S ll) of the ionized gas. Both Te(0 ill) and A''e(S ll) are 
systematically higher in low- mass galaxies with high EW(H/3). 
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Fig. 14. — (a) Oxygen abundance 12 + log 0/H vs. stellar mass M^, for the LCG sample. 
The symbols are the same as in as in Fig. |2l The linear least-squares best fit to the data 
is shown by the solid lin e. The fits to the green pea sample and star-forming SDSS galaxies 
of lAmorm et al. I ( l2010l ) are shown respectively by dashed and dash-dotted lines, (b) The 
N/0 vs. M^, relation for the LCG sample. The symbols are the same as in Fig. [2j 
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Fig. 15. — (a) The star formation rate SFR(Ha) and (b) the specific star formation rate 
SSFR(Ha) vs. the total stellar mass M*. The dashed lines in (b) correspond, from left to 
right, to loci with SFR = 1, 10 and 100 Mq yr~^, respectively. 
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Fig. 16. — Age t(old) of the oldest stellar population in the galaxy vs. its total stellar mass 
M*. The dashed lines from left to right correspond to loci with the SFR = 1, 10 and 100 
Mq yr~^, respectively. 



